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Abstract 

We calculate the fidelity of transmission of a single qubit between 
distant sites on semi-infinite and finite chains of spins coupled via the 
magnetic dipole interaction. We show that such systems often perform 
better than their Heisenberg nearest-neighbour coupled counterparts, 
and that fidelities closely approaching unity can be attained between 
the ends of finite chains without any special engineering of the sys- 
tem, although state transfer becomes slow in long chains. We discuss 
possible optimization methods, and find that, for any length, the best 
compromise between the quality and the speed of the communication 
is obtained in a nearly uniform chain of 4 spins. 

1 Introduction 

Quantum Information Processing (QIP) offers several advantages over classi- 
cal computation in solving problems associated with large and/or dynamical 
systems. One of the main goals of current research in the field is to find 
ways and means of reliably transmitting quantum information, encoded in 
quantum bits (or qubits), over arbitrarily large distances. In order to en- 
sure the information is received as the sender intended, the qubits must be 
protected from interacting with the environment in any way. For this rea- 
son, until now the qubits of choice have mostly been photons, which have 
an extremely small interaction cross section and can thus be used to gener- 
ate quantum states that are sufficiently robust to perform protocols such as 
quantum cryptography and teleportation ( P, Ej [3l HI [3, [6] ) . 
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However, in recent years much effort has been dedicated to studying sys- 
tems in which quantum information is encoded in stationary qubits, and is 
propagated from one part of the system to another by the interaction between 
the system's components. One of the simplest geometries in which this can be 
achieved is a one-dimensional chain of interacting particles, where the qubit 
is encoded in some internal degree of freedom (which we call 'spin', using |0) 
for spin 'up' and |1) for spin 'down'). Originally an exchange-coupled chain 
of spins with a constant nearest-neighbour (nn) interaction was studied [7] ; 
it was found that a qubit could be transferred with a fidelity exceeding the 
maximum classical value in a time that grows polynomially with the length 
of the chain. Subsequently, it was shown that such simple chains allow trans- 
mission fidelities arbitrarily close to unity also if the qubit is taken to be a 
carefully designed 'wave packet', provided the sending and receiving parties 
can access a sufficiently large portion of the chain [8]. Stronger results can be 
found for more complex systems: in the absence of structural imperfections, 
an XY Hamiltonian on a hyper-cubic lattice allows perfect state transfer 
[HI [ini [H] , as does a pair of parallel spin chains [12] , or a spin chain acting 
as a quantum wire connecting two qubits [13] . Plenio et al. [T^ have studied 
the situation for chains of harmonic oscillators (i.e. where each particle on 
the lattice possesses a continuous, rather than a discrete, degree of freedom), 
while Hartmann et al. [I5] have recently found that quantum information 
can be made to propagate with arbitrarily high fidelity through both oscil- 
lator and spin chains near a quantum phase transition, provided the ground 
state and the lowest excited state of the system are not degenerate. However, 
this transfer is exponentially slow; more rapid transmission is possible at the 
quantum critical point, but at some cost to the fidelity. 

It has also been shown that high fidelities can be attained by engineering 
the strength and the nature of the interactions between the spins [TB]. How- 
ever, this requires structures that would be very difficult to manufacture, 
both because the communicating parties would need to have an extremely 
high degree of control over the system, and because the component spins 
would have to exhibit nn couplings only, with precisely defined strengths. 
This is clearly an idealization, because long-range interactions are also likely 
to be present. Previous work has been dedicated to systems more 'realistic' 
from this point of view. Kay et al. have studied finite spin chains in which 
the total Hamiltonian accounts for the presence of local magnetic fields and 
a potential of the form of the magnetic dipole interaction [17]. The approach 
adopted in this case was to pre-determine a spectrum of eigenvalues that 
would ensure perfect state transfer for all chain lengths, and subsequently 
derive the corresponding local fields and the inter-spin distances by solving 
an inverse eigenvalue problem. This method is applicable to a Hamiltonian 
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containing any number of parameters, and could provide very useful theo- 
retical guidelines on the optimal way to structure a system, although, once 
again, the gap between theory and experiment may prove difficult to bridge. 

In this work, by contrast, we investigate simple one-dimensional arrays 
of spins interacting via a pure magnetic dipole interaction. We allow no 
site-specific locally-tunable fields; nevertheless, we show that fidelities for 
quantum state transfer closely approaching unity can be attained between 
the ends of finite chains, without any special engineering of the system. Fur- 
thermore, because of the long range interaction, the transfer rate grows poly- 
nomially in the system size, rather than exponentially, as in the case studied 
by Hartmann et al. [15]. Our results may be relevant to two- level atoms 
in atomic traps [IHl [IH] or to one-dimensional arrays of endohedral fuUerene 
species encapsulated within carbon nanotubes [20] , as well as to natural mag- 
netic dipolar systems such as LiHoF4 [21], and finite spin chains of a more 
complex design, for example engineered from arrays of quantum dots |22j . 

2 The System 

We build on the work done by Bose [7] on transferring quantum information 
through an infinite, uniform chain exhibiting isotropic nn interactions only. 
In this system, the qubit is represented by a single fiipped spin, which prop- 
agates between different sites in a manner defined by a time-independent 
Hamiltonian of the form: 

J ^ 
H = -^S,^i,Y.^'-^'-BY.crl (1) 

where N is the number of spins in the chain, a' = (o"^, cj,, o"^, ) are the 
Pauli spin matrices for the z*'' spin, i? > is a uniform magnetic field, and 
(J/2)5j_i_ij is the coupling strength between spins i and j, which is non- 
zero for nearest-neighbouring spins only. If the ground state of the system is 
expressed as | i) (8) I i) ••• I i), it has been shown that, provided the chain 
is sufficiently short, perfect or near perfect state transfer can be achieved 
by simply letting an initial state of the form | t) (8> I i)(8>---<8)| i) evolve 
naturally in time according to the effects of H. It is important to note that 
[H, J2iLo (^l] = 0, that is, the Hamiltonian conserves the total magnetization 
M of the system, allowing the chosen initial state to evolve only into states 
in which one spin is flipped at any given time. 

We propose to investigate the quality and efficiency of quantum state 
transfer through infinite and finite chains of spin-| fermions coupled by long- 
range interactions having the form of the magnetic dipole interaction. These 



3 



two systems differ in that the finite chain has end points, whereas the infinite 
chain does not; indeed, due to the periodic nature of the infinite chain we 
will hereafter refer to it as a ring. We examine a simplified system in which 
any external magnetic field is constant and parallel to the axis joining the 
dipoles, which is chosen to coincide with the z direction. The alignment of 
the magnetic and dipole axes is not a trivial point, but a necessary condition 
to ensure that the total magnetization M remains a good quantum number, 
and allows us to work in the sub-space where only one spin is flipped with 
respect to the ground state, reducing the Hamiltonian from a 2^ x 2^ to a 
NxN matrix. Within this subspace the effect of the magnetic field is to add a 
constant to the energies. We will hereafter omit this constant. Following the 
notation used in [7], we denote with |0000...0) the (unique) ground state of 
the system (i.e. all spins facing down, parallel to the external field) and with 
\j) the block of states in which the spin at the j^'^ site has been flipped from 
to 1. For simplicity, we assume there are no thermally excited spin-flips in 
the system. We adopt a Hamiltonian of the form: 

if, = ^[Si.Sj-3SrS|] (2) 

where C is a constant. Si and Sj are the total spin operators at sites i and j, 
and Sf and S? are the respective z components. The value of C is determined 
by the type of particle in the chain. For a system of spin-^ fermions (e.g. 
electrons) C is given by: 

where /iq is the permeability of free space, /x^ is the Bohr magneton, g is 
the electronic Lande g-factor and h = h/2'K. Throughout this paper we will 
assume that /io = /iB = ^ = l, so that: 

We define a to be the spacing between neighbouring fermions. In this case, 
the strength of the interaction between nearest neighbours is: 

(z|if,K±l) = ^ (5) 

For the results shown we define our length, energy and time units by setting 
the nearest-neighbour separation and the interaction energy between nearest 
neighbours to unity. However, eqn.([5]) implies the Hamiltonian has an overall 
scaling factor of so a uniform compression or expansion of the system 
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should only have the quantitative effect of re-scaling the system's energy by 
a constant. Therefore, provided the chain remains uniform and the number 
of component spins is fixed, the energy and performance of a chain of any 
size can be extrapolated by simply adjusting the value of a as necessary. 



2.1 Rings 



We initially consider a ring of N spins in its ground state. Our aim is to 
calculate the maximum fidelity of transmission of a qubit from site r to a 
distant site s, as a function of time and number of spins in the ring. We 
denote the initial and final states of the system by |r) and \s) respectively. 
The expression for the maximum fidelity of quantum state transfer is given 
by [7] as: 



3 6 2 

where f^^ (t) is the propagator, which is calculated from the following: 



(6) 



N 



■iEmt 



(7) 



m=l 



We assume that the eigenvectors |m) of the system can be expanded using 
the basis formed by the \j) states. Imposing periodic boundary conditions 
allows us to express these eigenstates as Bloch states, so that for a ring of 
circumference L = Na: 



\m) 



N-l 



(8) 



The propagator can then be re-expressed as: 



N 



&) = ^ E ^ 

m=l 



ik,na{r-s) ^-iEmt 



1 N-l 
^ m=0 



Using eqn. ([8]), we calculate the energies of the system, which can be written 
as: 



1 N-l 

Em = {m\H,\m) = ^ E e^'^^^^'-^^^l^^dlj) 



(10) 



We require only matrix elements for which i ^ j, since by symmetry the 
diagonal terms of are independent of TV, and therefore change the energy 
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of the system by a constant shift. For evenly spaced spins with nearest- 
neighbour separation a, these off-diagonal terms are: 

As eqn. ( ITTI) only depends on the difference |j — i|, we can reduce the double 
summation to a single sum over j by fixing a value of i. For convenience we 
choose i = 0, so that: 



N_ N_ 

F -i- V r^^wa 1 _ C ^ COS {kmja) , . 




Figure 1: (Colour online) Figure (a) shows the maximum fidelity that can be achieved in 
transferring an input state |1) to an output state j-y + 1) or | ^2^ ) on a ring of spins coupled 
by dipole-dipole (dashed blue curve) or nearest-neighbour (black curve) interactions, as a 
function of N. The red dotted line at F = 2/3 indicates the highest fidelity for classical 
transmission of a quantum state. Figure (b) shows the time at which the fidelity first 
peaks in this system. Units are as specified in Section [2] 



Fig. 1 (a) shows the maximum fidelity of state transfer for rings of = 3 



to = 30 spins, when the sending and receiving parties are located at 
diametrically opposite sites, for which |r) = |1) and |s) = |y + 1) or \s) = 
' ^+-) for even and odd N respectively. We have assumed that transfer occurs 
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along the arc joining sites r and s, though this may not be exactly the case. 



Fig. 1 (a) also shows the performance of a ring in which the spins are coupled 
by nn interactions only. We note that for A^ > 3 the performance of the 
nn-coupled ring is slightly better unless A^ = 6 or A^ = 12. We also find 
that the times of optimum transfer tend to rise as we increase N in both the 
dipole and the nn-coupled rings. 
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2.2 Single Qubit Transfer in Uniform Chains 

We now extend the previous analysis to a finite chain, calculating the full 
Hamiltonian of the system, which has the form: 



c 



c 



E 



2a3|j -i|3 
1 



(13) 
(14) 



where (000..0|ifrf|000..0) = — ^ S(fc,o Jkh^ ground state energy of the 

system. 



The fidelity of state transfer between the ends of the chain is obtained 
by taking |r) = |1) and \s) = \N). We find that F/^(t) exhibits three 
'trademark' features. 




5000 10000 15000 20000 



Figure 2: The evolution in time (abseissa) of the fidehty of state transfer between sites 1 
and 10 of a uniform chain of magnetic dipole-coupled spins. We note the regularity of the 
oscillation and the high value of F{raax). Units are as specified in Section [21 



First of all, the maximum value of F{t)i^ is close to unity. Secondly, the 
value oiFij^{t) oscillates between 1/2 and the maximum (which we call Fmax) 
with a regular frequency, which is generally quite small, implying that state 
transfer occurs slowly (fig. Finally, the period of oscillation of Ff^{t), 
which we call T, is uniquely defined by the energy splitting AA between the 
two lowest eigenvalues of Hd{N). The transfer process is therefore dominated 
by the beating of two nearly degenerate states localized near the ends of the 
chain. This behaviour is explained by the variation of the on-site energies of 
the spins as a function of j, shown in fig. [HI it is immediately evident that the 
most favourable positions for a spin to fiip are sites 1 and A^. Consequently, 
states |1) and |A^) are the most strongly coupled to the system's (two) bound 
states, which are shown in fig. [H In this system, this phenomenon is a natural 
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consequence of the geometry, but systems in which the spin flip energy is 
specifically chosen site by site have also been studied [23] ■ 



<j |H| j> 
-llr 

-11.5- 

-12 - 
-12.5- 

-13 - 
-13.5- 
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Figure 3: The on-site energy as a function of the site j of the spin flip for a chain of 15 
magnetic dipole-coupled spins. We note that the energies at sites 1 and 15 are much lower 
than the rest. Units arc as specified in Section [21 



litt 




Figure 4: (Colour online) The energy splitting AEm between the ground state |0000..0) 
and states with a single flipped spin for a uniform chain of N spins, showing the evolution 
of the two bound states. Values of AEm of the same index m, counting from the bottom 
of the spectrum, are shown in the same colour. Units are as specified in Section [51 



The period of F{^(t) is related to AA by: 



T 



2n 
AA 



(15) 



Consequently, the time at which Ff^{t) first peaks is: 

T n 

t[Fmax) =2 = AX 

This time rises with chain length, as the splitting AA decreases with increas- 
ing N. 



(16) 
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A summary of our results for = 2 to = 23 spins is shown in Fig. 5(a) 
We note that, in addition to A^ = 2, A^ = 3 and A^ = 4 also give perfect 
transfer, and in general Fmax > 0.9. This is a marked improvement on the 



performance of a nn-coupled chain (also shown in Fig. 5(a)); in particular, it 
seems that by replacing the nn couplings with dipole couplings we no longer 
obtain poor transfer when A'^ is a multiple of 3 [7]. Unfortunately, we again 
observe a 'trade-off' between fidelity and time, which becomes particularly 
evident for A^ > 6. In fact, we find that at large A^, t{Fmax) goes as the cube 
of the chain length (fig. 5(b)). It is therefore evident that in long chains it 



will take an impractical length of time to complete the protocol unless the 
system can be optimized in some way. 



Logiot 



15 20 





Figure 5: (Colour online) Figure (a) shows the maximum fidelity that can be achieved 
in transferring an input state |1) to an output state |A^) in a chain of spins coupled by 
dipole-dipole (dashed blue curve) or nearest-neighbour (black curve) interactions, as a 
function of N. The red dotted line at = 2/3 indicates the highest fidelity for classical 
transmission of a quantum state. We note that the dipole-coupled chain almost always 
performs better. Figure (b) shows the time at which the fidelity first peaks in these two 
systems, plotted on a Logio scale. The red dotted curve is the function y — . We see 
that at large L the dashed blue and red dotted curves are parallel, indicating that the 
transfer time scales as the cube of the chain length. Units are as specified in Section [2l 



2.2.1 Structural Optimization 

We now analyze the efficiency of state transfer for a fixed chain length, as a 
function of the number of spins in the chain. We define r as the transmission 
time giving maximum fidelity at unit chain length, i.e.: 

- = (17) 

Fig. [6] shows a plot of function of N, which reveals two interesting 

features. The first is the presence of a minimum at A^ = 4, which we will 
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discuss subsequently. The second is the fact that r tends to a constant 
at large L. This indicates that, above a certain threshold value of N, the 
evolution of the system is determined almost exclusively by the magnetic 
dipole coupling between spins 1 to g with spins N — q to N, irrespective 
of the number of spins that separate these two 'clumps'. To explore this 
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Figure 6: The behaviour of r as a function of the number of spins in the chain. Note the 
minimum at iV = 4 and the flatness of the curve for N > 15. Units are as specified in 
Section H 

hypothesis, and determine the behaviour of r for large N, we work with 
states \B) and \E) localized at the beginning and end of the chain, which 
are the bound-state eigenfunctions of a semi-infinite chain extending to the 
right and the left, respectively. We take: 



n=l 



B) = J2(^n\n) (18) 



and 



\E) = Y.an\N+l-n) (19) 



n=l 



The energy splitting of the two lowest eigenvalues of HalN) in a finite chain 
is: 

AA = 2{B\Hd\E) (20) 

From dn]): 

{t\H,\j) = H,i\t-j\) (21) 

Hence: 

{B\Hd\E)= ^ a*^am{n\Hd\N + 1 - m) = ^ ala^Hd{\N + 1 - m - n\) 

n,m=l n,m=l 

(22) 
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We adopt a 'dummy' variable X = |z — j|, so that: 



2a3X3 



dX 2a^X^ 

Using ([22]) and the fact that L = a{N — l), we can expand Ha{\N + l — m 
as a Taylor series to first order in 6 = m + n — 2. Then: 



Hence: 



with: 



dX 



C ^ 3Ca{m + n - 2) 



X=L 



2L3 



2L4 



(23) 

(24) 
-n\) 

(25) 



2" 

1 
~2 



1 a 

— a*^am + jj 3<a„(m + n-2) 

n,m=l n,m=l 

Q aR 



Q= E 



(26) 
(27) 

(28) 



n,m=l 



R= E 3a*am(?7i + n — 2) 

n,m=l 



(29) 



lOt 
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Figure 7: (Colour Online) Comparison between the predictions of our q-spin model (red 
dashed points) and the data calculated from the treatment of the system in its entirety. 
We note the model becomes increasingly accurate at large N. 

We find it is possible to model the asymptotic behaviour of the system 
very accurately if we assign to the coefficients a„ the amplitudes of the ground 
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state eigenvector of the q x q sub-matrix generated by truncating the full 
Hamiltonian Hd{N) for an arbitrarily large N (fig. [7]). The values of a„ used 
to obtain the fit in fig. [7] come from the ground state eigenvector of the 4x4 
sub-matrix of ifrf(14), and give Q ~ 0.325 and R ~ —0.957. These quantities 
show only a very weak dependence on A^, so we have assumed them to be 
constants. The fact that Q < 1 shows that the transfer rate for chains of 
many spins is always less than that attained between two completely isolated 
spins; equation ( |28l) indicates this is a result of interference between positive 
and negative components in the localized states \B) and \E). 

Therefore, only chains with few spins can improve on the performance 
of a simple dipole pair, as shown by the minimum in the function t{N) 



at N = A (fig. [H]). This result, together with fig. 5(a), shows that, in a 
uniform chain, the best compromise between the quality and the speed of the 
communication is obtained with 4 spins. This occurs because for short chains 
the bound states at the two ends have a large overlap, i.e. there exist terms 
in eqn. (!22|) which simultaneously have significant positive values of a* and 
small values of + 1 — m — n|. We have attempted to optimize the uniform 
4-spin chain still further, and find it is possible to improve its performance 
slightly by modifying the positions of the inner spins while maintaining mirror 
symmetry. For a chain of unit length, this corresponds to taking ri.2 = ^3^4 ^ 
0.314 and r2.3 ~ 0.373, which yield a value r ^ 0.512. However, a comparison 
with r ~ 0.568 for a uniform chain (ri 2 = 'r2_3 = r3 4 ^ 0.333) shows the 
improvement is minimal. 



2.2.2 Input and Output Optimization 

We now investigate the effects of altering the initial and final states |r) and 
|s), while leaving the structure of the chain intact. 

If the starting and ending points are chosen at random, the characteristic 
oscillation of F^^(t) is lost, unless either |r) = |2) and \s) = \N), or |r) = |1) 
and \s) = \N — 1). However, in both cases the signal is considerably noisier, 
and the maximum fidelity is greatly reduced. This is a result of the lesser 
efficiency of coupling to the bound states as one moves away from the ends 
of the chain (cfr. fig. H]). 

Conversely, it is possible to boost the maximum fidelity to unity and 
smooth out all noise in the signal by encoding the states |r) and |s) in two 
or more adjacent spin^(fig. [H]) . This is equivalent to adopting |r) = \B) and 
\s) = \E), where the a„ are now obtained from the first and last n coefficients 
of the ground state eigenvector of Hd{N), with the additional condition that 

""^The possibility of improving state transfer by encoding a state in more than a single 
spin is also discussed in [S]. 
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Figure 8: The evolution in time (abscissa) of the fidehty of transmission of an input state 
of the form Ci|l) +C2I2) to an output state of the form Cg\9) + Cio|10) in a uniform chain 
of 10 dipole-dipolc coupled spins. Comparing with fig.[H note that the curve is smoother 
and the maximum fidelity has increased. Units are as specified in Section [51 

J2n knP = 1- This choice of input and output states leaves the transfer time 
unaffected. 

3 Discussion and Conclusions 

We present a scheme for transferring quantum information through infinite 
and finite chains of spins coupled via a pure magnetic dipole interaction. This 
differs from much previous work in that the dipole interaction is long-range, 
making for a system in which every spin interacts with all other spins in the 
system, rather than with nearest neighbours only. We find that, in general, 
the maximum fidelity achievable by using a dipole-coupled system to transfer 
a state between two maximally distant sites is greater or equal to that which 
can be attained in a system exhibiting nearest-neighbour interactions only. 
The finite chain, in particular, can be engineered to give unit fidelity by 
simply adjusting the placement of the spins and the input and output states. 
We have verified this result only for L = 2 to L = 23, but believe it extends 
to longer chains also. 

The main weakness of both our systems is length of time taken to com- 
plete the protocol, which increases polynomially in the size of the system. 
However, we find that for a finite chain this obstacle can be considerably 
lessened by simply modifying the relative placement of the spins. Therefore, 
it does not necessarily preclude the possibility of being able to transmit infor- 
mation over longer distances on useful timescales. Furthermore, the protocol 
seems to be reasonably robust against errors in spin placement; if we define 
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a "failure rate" as the probability that the fidelity at time t{Fmax) will fall 
below the classical value, we find that, in a uniform 4-spin chain of arbitrary 
length, a random error of 2 % on the placement of each spin yields a fail- 
ure rate of approximately 5 %. The corresponding uncertainty on t{Fmax) is 
significantly greater, but as F^^{t) is a slowly varying function of time, it is 
quite unlikely that the fidelity at t{Fmax) =t ^ti^Fmax) will have fallen signifi- 
cantly below the maximum. Therefore, the simple and predictable behaviour 
of the fidelity in time in a finite chain greatly increases the probability of car- 
rying out successful state transfer. The long range interactions also open up 
the realistic possibility of measurements on individual spins, and it would be 
interesting to investigate in the future whether this can increase the speed of 
quantum state transfer as in Ref. 
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